A practical formula for the radiated angular momentum 
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We present a simple formula for the radiated angular momentum based on a spin- weighted spheri- 
cal harmonic decomposition of the Weyl scalar ifti representing outgoing radiation in the Kinnersley 
tetrad. We test our formula by measuring the radiated angular momentum from three simulations 
of non-spinning equal-mass black-hole binaries with orbital angular momentum aligned along the a;, 
y, and z axes respectively. We find that the radiated angular momentum agrees with the differences 
in the remnant horizon spins and the initial angular momentum for each system. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.Bw, 04.30.-w 



I. INTRODUCTION 



Recent breakthroughs in numerical relativity [j], 0, HJ 
have made it possible to accurately evolve black-hole — 
black-hole 0, [f| d, 0] and black-hole — neutron-star d, [§] 
binaries for many orbits. The feasibility of evolving 
generic black-hole binaries was demonstrated in [Id] 
where unexpectedly large recoil velocities of the merger 
remnant [lfjLflll] were found. One of the important meth- 
ods to measure the accuracy of the full numerical evolu- 
tions is to monitor the conservation of the energy, linear 
momentum, and angular momentum. A practical way 
of doing this in black-hole binary systems is to compute 
the final black hole mass, linear momentum, and angular 
momentum from the (isolated) horizon expressions, as 
well the radiated energy, linear momentum, and angular 
momentum as measured far away from the system (for- 
mally at null infinity, ,y + ). The sum of the corresponding 
remnant and radiated quantities should converge to the 
initial mass, linear momentum, and angular momentum 
of the system. 



Expressions for the energy and momenta of an isolated 
system have long been established (see Refs. [12J, ll3j for 
reviews) . Currently, the most popular method to extract 
gravitational radiation is to compute the Weyl scalar if>± 
in a quasi-Kinnersley frame or a frame that reduces 
to one asymptotically. For practical purposes, it is then 
convenient to express the radiation formulae for the mass, 
linear momentum, and angular momentum in terms of 
-04. That was done in Ref. [l5[, see formulae (22)-(24) 
there, where the radiated angular momentum was only 
given along the z-axis. We extend the computation of the 
radiated angular momentum to all the three Cartesian 
components. These are useful to verify the conservation 
of the angular momentum in generic black- hole binaries, 
for instance, when spin precession is present [lol Ull. fl^|. 



II. DERIVATION 

The angular momentum on a space-like slice S + of J !+ 
(e.g. the intersection of a u = t—r = const slice with J^ + ) 
is given by [l3[ 

J[i] = Re U dS^q A [2V>? - 2a°<5a° - B(a°a )] 

(1) 

where the superscript indicates the coefficient of the 
leading-order term in an expansion in 1/r, the spin- 
coefficient a and Weyl scalars ijiQ, • • • } "04 are given in 
terms of a null tetrad (Z° , n a , m a , m a ) by (Eqs. (3.7-3.8) 
of Ref. [13) 
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(C abcd l a n b l c n d - C abc dl a n b m c m d ) , (5) 
-C abc dn a fh b n c m d , (7) 
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(where C a bcd is the Weyl tensor), and have the asymp- 
totic form (Eq. (4.4) in 



17 


= a°r- 2 ^ 




(8) 


ipo 


= i>° r- 5 - 




(9) 


Ipl 


= ^r" 4 - 


K0(r- 5 ), 


(10) 


1p2 






(11) 




= ^° 3 r- 2 - 


(-0(r- 3 ), 


(12) 


i>4 


= y^lr- 1 - 


KO(r- 2 ), 


(13) 



the coordinates on the slice S + are such that the vol- 
ume element has the standard form for a 2-sphere, dS = 
sin6d9d(j), the capital A index ranges from 1 to 2 over 
the angular coordinates in the slice, the dyad q A is the 
standard unit-sphere dyad (q A — (1, z/sin#) in (9, (ft) co- 
ordinates) , the index i ranges from 1 to 3 and represents 
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the three Cartesian components of the angular momen- 
tum, an overbar indicates complex conjugation, the three 
£^ are rotational Killing vectors intrinsic to E + , and the 
3 and 5 operators acting a function v s of spin-weight s 
have the form 



-d e v s 



sin 9 



d^Vg + sv s cot 9, 



Qv s = -d e v s + -^—rd^Vs - sv s cot 9. (14) 
sm9 

The tetrad used in Eq. (TT]) is partially fixed by the 
choices l a = V a u, u = (t — r)/\/2 + O(logr), and 
m a = (0,0, l,i/sin0)/\/2 + C>(l/r), where r is a lu- 
minosity distance (we consider here the case where the 
spacetime and coordinates are such that the metric ap- 
proaches diag(l,l,l) at large r and hence the coordinate 
r is a luminosity distance). Note that l a must therefore 
be hypersurface orthogonal. The quasi-Kinnersley l a , by 
contrast, is not hypersurface orthogonal when the spin 
of the hole is non-vanishing and will give an incorrect 
value for the angular momentum (e.g. zero for a Kerr 
black hole). However, as we will derive the radiated an- 
gular momentum in terms of 1(14, and the Kinnersley (l8j 
l a differs from V a it by a type II rotation (which leaves 
-04 invariant) and a trivial boost of V2, the Kinnersley 
i/>4 can be used for the radiated angular momentum for- 
mulae derived below since it is invariant under type II 
rotations. 

Note that for any spin-zero function / and Killing vec- 
tor (of the unit-sphere metric) £ , 



Re 
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dS£ A q A dfj = Re I j>^ + dSeU ^(f - /)] 

(15) 

i.e. only the imaginary part of / contributes to the inte- 
gral. Thus the final term in Eq. ([1]) does not contribute 
to J. 

Let us define = £j4 q A where £4 are the standard 
rotational Killing vectors about the x, y, and z axes re- 
spectively, then 



= icos0cos^-sin0 = 2i A /-(_ 1 Y 11 -_ 1 Y 1 _ 1 ), 



= icos#sin0 + cos0 = 2\\—{-\Y\ \ + -\f\ _i), 



<I> -isin6» = -2i\j — -iY 10 , 



(16) 



where the spin-weighted spherical harmonics have the 
form [2l| 



sYtm = (-iyJ |^yy 3 s Y lm (for a > 0), 



sY tm = Jj^yy s " s Y tm (for s < 0). (17) 



Note that all three $rji can be expressed as a sum of 
£ = 1 spin-weighted spherical harmonics (of spin-weight 
-1). Hence the angular momentum loss is given by 

j w = -rj: Re jy dS$ [fl (2-0? - 2ct 3ct° - 2ct 3ct ) ) 

(18) 

where a dot indicates differentiation with respect to u. 
Using the following relations [13], 

tpl = -30>° + 2ct°3ct ) , (19) 

^ 2 °-0 2 ° = BV - 3 2 a° + <r°<r° - (20) 

< = (21) 
Eq. (fl"8|) reduces to (where we used result (fT5|) ) 

= ^ Re jj^dS^S^-BBV- 
SaW + 3a°3^° - a°3d° + S°3a )} . (22) 

Note that the first two terms do not contribute to Jm 
since 3>[j] contains only £ = 1 modes and cr° (which has 
spin-weight 2) cannot contain £ = 1 modes. Hence we 
have 

j M = Y^Rej^" d5$ H (-3CT 3CT° + 3CT 3CT°- 



a 3d +S°3a )} 



(23) 



In terms of the more standard choice of tetrad with l a — 
V a (t— r), the radiated angular momentum takes the form 

= — Re II dS$ H] (-3Nm + 3HBN- 
dt 16-k Lje+ 

HW + NdH)}, (24) 

where dtt-ff = dtN = -04- 

After expressing H and AT in terms of spin-2 spherical 
harmonics 



H — ^ Hi m -lYlm, 
l,m 

N = ^N lm2 Y lmi 



(25) 
(26) 



Eq. (|24|) becomes 

^ = "1^ E M^m^ m ,41,, m ,}, (27) 
where 



a iml>m> ~ C lml'm> ~ C t' m' £rn^ 



(28) 



and 



dS*w$y/(l + 2){t-\) xYtm iY i>ml 
■y/(£' -2)(t' + 3) 2 Y em 3 Y t , m ,). (29) 
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The coefficients d 



A 



have the simple form 



2iy/£(£ + 1) - ram' 5| m _ m /|i, 



2(m - m')y/l{i + 1) - mTO'5| m _ m /| 1; 



(30) 



[22| . The radiated angular momenta then has the form 



f m 



± v /^±l)-ra(m-l) Im{ir /fll JVim-i}J , (31) 
^- V (\/^±l)-m(ra±l) Re{i2, m iV> m+ i} 

Im 

-y/t{t + 1) - ra(ra - 1) Re{H/ m iV< m _i}) , (32) 



(33) 



Note that if -04 



2^, 



then Ne m = 

(-ir/Ioo^-mWdr, and Jfc m = Jl^NtMdT. In 
practice one needs to integrate from the starting time 
of the simulation rather than t = — oo. This presents 
no significant problem for black-hole binary simulations 
as the error introduced by assuming no radiation in the 
past is of the same order as the error introduced by us- 
ing initial data without the correct radiation content. In 
addition, these simulation typically start when the ra- 
diation is quite small compared to the late-inspiral and 
merger waveforms. 



III. NUMERICAL TESTS 

We tested the newly introduced formulae for the x, 
y, and z components of the radiated angular momen- 
tum by evolving the non-spinning, quasi-circular SO con- 
figuration of Ref. [l!| rotated such that orbital angu- 
lar momentum is aligned along the x (SOX), y (SOY), 
and z (SOZ) axes respectively. For reference we pro- 
vide the initial data parameters in Table U for these 
configurations. The total radiated angular momentum 
for each configuration is reported in Table HI We ran 
the SOX and SOY configurations with our new AMR 
code [l(| with 9 levels of refinement and central res- 
olution of M/40. We extracted the waveform at radii 
r = 25M, 30M, 35M, and 40M. We fit the calculated 
radiated angular momentum to a linear and quadratic 
polynomial in 1/r. We take the quadratic extrapolation 
to be the value at r — oo and the quoted errors are the 
differences between the linear and quadratic extrapola- 
tions. Note that we did not re-run the SOZ configura- 
tion, but rather re-analyzed the waveform from an exist- 
ing unigrid run with a relatively low central resolution of 
h = M/22.5 and extraction radii at r = 20M, 25M, and 



TABLE I: Initial data and total radiated angular momentum 
for quasi-circular, equal-mass black-hole binaries. The total 
angular momentum is J/M 2 = 0.8764 and the proper hori- 
zon separation is l/M — 10.01. The punctures are located 
at ±{X,Y, Z), with mass parameter m p /M — 0.4848, and 
momentum ±(P x ,P y ,P z ). The orbital parameter (X,Y,Z) 
and (P x ,P y , P z ) are obtained by the 3PN equations of motion 
for equal-mass binaries with orbital frequency Mfl — 0.0500. 
The puncture masses are obtained by setting the ADM mass 
of the system to (1.0000 ± 0.0005)M. 



Config 


SOX 


SOY 


SOZ 


X/M 


0.000 


0.000 


3.280 


Y/M 


0.000 


0.000 


0.000 


Z/M 


3.280 


3.280 


0.000 


P x /M 


0.0000 


0.1336 


0.0000 


Py/M 


-0.1336 


0.0000 


0.1336 


P z /M 


0.0000 


0.0000 


0.0000 




0.236 ± 0.007 


< 10" 7 





< 10" 6 


0.236 ± 0.007 





< 10" 6 


< 10 -7 


0.24 ±0.01 



30M. In all cases we summed over all modes with i < 4 
[23| . The remnant hole for these configurations (SOX, 
SOY, and SOZ all represent the same binary) has spin 
5 remnant /Af 2 = 0.639 ±0.001 indicating that the angular 
momentum loss is S J/M 2 = 0.237 ± 0.001, which agrees, 
to within the error estimates, with the radiated angular 
momenta reported in Table HI 



IV. CONCLUSION 

We used the Linkage definition of the angular momen- 
tum to derive a formula for the radiated angular mo- 
mentum (|24|) . (I14p . (fT6f in terms of ip4. The remaining 
radiation formulae for the energy and linear momentum 
have already been given in this form in Ref. [H, Eqs. 
(22)-(23). 

We then derived a simple formula for the radiated an- 
gular momentum based on a spin- weighted spherical har- 
monic decomposition of ip±, and have demonstrated that 
the formula provides the correct radiated angular mo- 
mentum for binaries that radiate along all three Carte- 
sian directions. We also compared those results with the 
spin of the horizon of the final black hole minus the total 
ADM angular momentum of the system, and find very 
good agreement. Note that we also used these formu- 
lae to produce this kind of check in a generic black-hole 
binary simulation involving unequal mass and unequal 
(precessing) spins, the SP6 run in Ref. and obtained 
satisfactorily agreement. This paper thus provides prac- 
tical formulae for direct application to extract radiation 
information from current numerical relativity simulations 
of compact sources. 

As a final point we note that there is an ambiguity in 
the definition of angular momentum for generic asymp- 
totically flat spacetimes [IH, [l3| that arises from the su- 
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pertranslation symmetry of This ambiguity can 

lead to inaccuracies in the computed angular momen- 
tum for generic gauges. In general one needs to confirm 
that the intrinsic metrics of the extraction spheres used 
to compute the modes of V>4 asymptotically approach 
ds 2 = r 2 d0 2 + r 2 sin 2 
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The definitions of and 3 in Eq. (U4[) differ from the 
conventions in [2(j by a factor of —1. We compensate for 
this factor in the definition of 3 Yi m in Eq. (JTTJ) . Thus the 
s Yg m used here are identical to those used in [20| . 
We have demonstrated that Eq. (|30|l is correct for 
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general. 

For these configuration the £ — 2 modes dominate the 
contribution to the radiated angular momentum. 



